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Abstract 

We show analytically that the dynamical equations for an underdamped 

ladder of coupled small Josephson junctions can be approximately reduced 

to the discrete sine-Gordon equation. As numerical confirmation, we solve 

the coupled Josephson equations for such a ladder in a magnetic field. We 

obtain discrete-sine-Gordon-like IV characteristics, including a flux flow and a 

"whirling" regime at low and high currents, and voltage steps which represent 

a lock-in between the vortex motion and linear "phasons", and which are 

quantitatively predicted by a simple formula. At sufficiently high anisotropy, 

the fluxons on the steps propagate ballistically. 

PACS numbers:74.40.+k, 74.50.+r, 03.20.+i, 05.45.+b, 85.25.C 
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The discrete sine-Gordon equation has been used by several groups to model so-called 
hybrid Josephson ladder arrays [[LfU . Such an array consists of a ladder of parallel Josephson 
junctions which are inductively coupled together, e. g. by superconducting wires |4]]. The 
sine-Gordon equation then describes the phase differences across the junctions. In an applied 
magnetic field, this equation predicts remarkably complex behavior, including flux flow 
resistance below a certain critical current, and a field-independent resistance above that 
current arising from so-called "whirling" modes f2|. In the flux flow regime, the fluxons in 
this ladder propagate as localized solitons, and the IV characteristics exhibit voltage plateaus 
arising from the locking of solitons to linear "spin wave" modes. At sufficiently large values 
of the anisotropy parameter r\j defined later, the solitons may propagate "ballistically" on 
the plateaus, i. e. may travel a considerable distance even after the driving current is turned 
off. 

In this Letter, we show that this behavior is all found in a model in which the ladder is 
treated as a network of coupled small junctions arranged along both the edges and the rungs 
of the ladder. This model is often used to treat two-dimensional Josephson networks || , and 
includes no inductive coupling between junctions, other than that produced by the other 
junctions. To confirm our numerical results, we derive a discrete sine-Gordon equation from 
our coupled-network model. Thus, these seemingly quite different models produce nearly 
identical behavior for ladders. By extension, they suggest that some properties of two- 
dimensional arrays might conceivably be treated by a similar simplification. In simulations 
0-0, underdamped arrays of this type show some similarities to ladder arrays, exhibiting 
the analogs of both voltage steps and whirling modes. 

Our ladder consists of coupled superconducting grains, the i th of which has order pa- 
rameter <&j = Q e l6i . Grains i and j are coupled by resistively-shunted Josephson junctions 
(RSJ's) with current iy, shunt resistance Rij and shunt capacitance CV,-, with periodic 
boundary conditions (see Fig. 1). 

The phases evolve according to the coupled RSJ equations Ti6i/{2e) = Vi, J2j MijVj = 
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If xt I Ic-J2j(R/ 'RijWi-Vj) -Hjihjl 'L) sin(%- Aij). Here the time unit is t = h/(2eRI c ), 
where R and I c are the shunt resistance and critical current across a junction in the x- 
direction (see Fig. 1); l? xt is the external current fed into the i th node; the spatial distances 
are given in units of the lattice spacing a, and the voltage V* in units of I C R. = 
—4ireCI c R 2 /h for % ^ j, and M u = — Yjj^i^Iij-, where C is the intergrain capacitance. 
Finally, = (27r/$ ) // A • dl, where A is the vector potential. We assume N plaquettes 
in the array, and postulate a current I uniformly injected into each node on the outer edge 
and extracted from each node on the inner edge of the ring. We also assume a uniform 
transverse magnetic field B = f(po/a 2 , and use the Landau gauge A = —Bxy. We solve 
these equations numerically using a fourth-order Runge-Kutta algorithm. 

We now show that this model reduces approximately to a discrete sine-Gordon equation 
for the phase differences. Label each grain by (x, y) where x/a — 0, . . . ,N — 1 and y /a — 0,1. 
Subtracting the equation of motion for 9(x, 1) from that for 9(x,2), and defining = 
\[0{x, 1) + 9(x,2)], x( x ) — [0(^)2) — Q(x,\)\, we obtain a differential equation for x( x ) 
which is second-order in time. This equation may be further simplified using the facts that 
A x , y -,x±i,y = in the Landau gauge, and that A x ^ x ^ = —A Xt 2- x ,i, and by defining the discrete 
Laplacian x( x + 1) ~ ^x( x ) + x( x ~ 1) = ^ 2 x( x )- Finally, using the boundary conditions, 
I ext (x, 2) = - I ext (x, 1) = I, and introducing (fi(x) — x( x ) — A Xt2 - x ,i, we obtain 

[1 - 77 C 2 V 2 ]/^ = i - [1 - rf r V 2 }^ - sin(v9) + 2^ x 

cos{^(x)-^(x + i)}sm{[(p(x)-(p(x + i)}/2}, (1) 

i=±l 

where we have defined a dimensionless current i = I/I cy , and anisotropy factors 2i] 2 = 
Ry/Rx, 2rf c = C x /C y , and 2r]j = I cx /I cy . 

We now neglect all combined space and time derivatives of order three or higher. Sim- 
ilarly, we set the cosine factor equal to unity(this is also checked numerically to be valid a 
posteriori) and linearize the sine factor in the last term, so that the final summation can 
be expressed simply as V 2 0. With these approximations, eq. (1) reduces to discrete driven 
sine- Gordon equation with dissipation: 
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f3(p + + sin(p) - r/jVV = i, (2) 

where (3 = AneI cy R^C y /h. 

To confirm the accuracy of this reduction, we have solved the coupled Josephson equa- 
tions on an 8 x 1 (N=8) ring. Fig. 1 shows the resulting IV characteristics with (3 = 33, 
r/j = 0.71, and several values of /. (V) denotes the space and time-averaged voltage dif- 
ferences across the y-junctions. There are two regimes. The first is the "flux flow regime" 
where (V)(I, /) is roughly proportional to / (up to about /=3/8) and to I. In this regime, 
for each /, (V)(I,f) also exhibits a series of voltage steps (see below). The second regime 
corresponds to "resistance steps" in which (V) = RI independent of /. This regime is 
dominated by the whirling modes, and is also discussed further below. 

Fig. 2 shows an expanded low-current regime for / = 1/8, (3 — 61, and several rj/s. 
The voltage steps are very prominent. The IV characteristics are hysteretic on each of the 
steps, as shown with broken lines for rjj = 1.25. Similar hysteretic steps are well known 
in numerical studies of the discrete sine-Gordon equation |]|§. The critical current for the 
onset of voltage varies from about 0.2J C at 77 j = 0.71 to ~ for rjj > 1, in the less discrete 
regime, similar to results obtained in ||. 

Soliton behavior. In the absence of damping and driving, the continuum version of eq. 
(2) has, among other solutions, the sine-Gordon soliton ||10|| , given by 



<p s (x,t) ~ 4 tan 1 exp |(x - v v t)/y/rfi - @v%, J 



(3) 



where v v is the velocity. The phase in this soliton rises from ~ to ~ 2-k in a width 

Fig. 3 shows the local phase difference tp(x,t) for the Josephson ladder at two currents 
in the flux-flow regime (I / I c = 0.22 and 0.41), as well as one in the "whirling" regime, 
I/I c = 0.8. (/ = 1/8, (3 = 61, rjj = 1.25). The first two show clear soliton-like behavior, 
namely, an increase of <p(x, t) in steps of ~ 2n over a time interval dk/v v , as predicted by the 
sine-Gordon equation. The passage of the kink is accompanied by ripples arising from phason 
excitations. Typically, an integer number of ripple periods is found between successive kink 



passages. The local voltage 0(0, t) (shown as broken lines) in the flux- flow regime is due to 
such sine-Gordon solitons, modified by coupling to phasons. The sharp peaks in the local 
voltage correspond to the vortex passage (we have confirmed this by snapshots of the vortex 
motion in our simulations), while the exponentially decaying smaller peaks correspond to 
the phasons which couple to the vortex. In the whirling regime at I/I c = 0.8, (V) oc /, and 
the local voltage oscillates sinusoidally in time. Each oscillation period again corresponds to 
the passage of a vortex through a given junction, as can be confirmed by the fact that the v v 
thus found is consistent with the independently computed (V), via the Josephson relation. 

The step positions in Fig. 2 are determined by a locking of the vortex motion to 
the phasons. The phason dispersion relation is determined by linearizing the left-hand 
side of eq. (2) with i = 0. The result is tp m (x,t) oc exp(— t/2(3)e l( - kmX ~ UJmt \ where 
ujm = ±v/T+ Ar/j sin 2 [k m j '2) /\fj3, as obtained previously by several groups [0,0. The al- 
lowed wave vectors k m are determined by periodic boundary conditions: k m = 2irm/N,m = 
0, 1, 2, 3,4,.. .. To obtain the locking condition, note that the vortex circulates the ladder 
with frequency uo v = 2irv v /(Na). A resonance will occur if there are an integer number of 
phason cycles per vortex cycle. This condition gives uj m = na^with n = 1,2,3,..., or 

1 = n - (4) 

won ,fi + Arfj sin 2 (nm /N) 

All the voltage steps we have found in Fig. 2 satisfy this condition. At rjj = 0.71, for 
example, we identify resonances in the range 3 < n < 15, m = 1, from a high-resolution 
IV characteristic and its derivatives. Presumably, the resonances corresponding to higher 
n are weaker because the phasons are damped, relaxing over a time 2(3. At larger rjj, we 
can identify only a few values of n. The values of n for each step can also be found by 
enumerating the number of phason wavelengths between successive vortex passages, as in 
Fig. 3. In the inset of Fig. 2 we compare the positions of the steps thus located to the 
predictions of eq. (4). In all cases, only the m — 1 mode is necessary to account for the 
resonances. 

The transition to the resistive state ("Region II") occurs at n min = 4,2,2,1 for rfj = 

5 



0.5,1.25,2.5,5. This can also be understood from the kink-phason resonance picture. To 
a phason mode, the passage of a kink of width df. will appear like the switching on of a 
step-like driving current over a time of order dk/v v . The kink will couple to the phasons 
only if dk/v v > tt/ux, the half-period of the phason, or equivalently 

i > VTT? = ™ (5) 

Vfov Vj Vj 
This condition agrees very well with our numerical observations, even though it was obtained 
by considering soliton solutions from the continuum sine-Gordon equation. 

The fact that the voltage in regime I is approximately linear in / can be qualitatively 
understood from the following argument. Suppose that ip for Nf fluxons can be approxi- 
mated as a sum of well-separated solitons, each moving with the same velocity and described 
by f(x,t) = YlfJifj, where <p>j = <p s {x — Xj,t). Since the solitons are well separated, we 
can use following properties: sin[X^ <Pj] = J2j sirup j and / (p>j(pidx oc 6^. By demanding that 
the energy dissipated by the damping of the moving soliton be balanced by that the driving 
current provides(oc / dxiip(x)), one can show that the Nf fluxons should move with the 
same velocity v as that for a single fluxon driven by the same current. In the "whirling" 
regime, the f- independence of the voltage can be understood from a somewhat different 
argument. Here, we assume a periodic solution of the form ip = JZ- fw(x — vt — j / f) 
moving with an unknown velocity v where <p w (£) describes a whirling solution containing 
one fluxon. Then using the property ip{x + m/f) = ip(x) + 2irm, one can show after some 
algebra that sin[X)j <p w {x — vt — j / f)} = sm[Nf(p w (x — vt)]. This means that Nfip w (x — vt) 
is a solution to eq. (2) as is (p w (x — vt). Finally, using the approximate property <p w (£) ~ £ 
of the whirling state, one finds v = v/(Nf), leading to an f-independent voltage. 

Ballistic soliton motion and soliton mass. A common feature of massive particles is their 
"ballistic motion," defined as inertial propagation after the driving force has been turned off. 
Such propagation has been reported experimentally [|11[] but as yet has not been observed 
numerically in either square or triangular lattices In the "flux-flow" regime at r]j = 

0.71, we also find no ballistic propagation, presumably because of the large pinning energies 
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produced by the periodic lattice. (The critical current for soliton depinning at rjj = 0.71 is 
about 0.2I C , about twice that calculated for a square lattice ||12| .) However, at r/j > 1, we 
do observe ballistic motion in the flux flow region I. As an example, Fig. |] shows V(t) and 
dV/dt at junction 0, for rjj — 5, I/I c — 0.41 (on the n = 1 voltage step), after the driving 
current is switched off at t = 0. The washboard-like ridges of V on a decreasing background, 
and perhaps more clearly, the spikes in dV/dt, correspond to a vortex passing through this 
junction. The increasing distance between peaks indicates that the vortex is slowing down. 
The vortex circulates at least five times around the ring before stopping - a fact which is 
also verified by direct observation of the vortex motion in real time. Qualitatively similar 
behavior was also observed for slower vortices on the lower current steps so long as rjj > 1. 
This behavior can be understood by noting that increasing t]j increases the width of the 
kink dk, and thereby effectively makes the ladder seem less discrete. The fluxon at large t]j 
therefore has a much lower depinning current and can propagate ballistically. This suggests 
the interesting experimental possibility that one can tune the vortex pinning current and 
effective mass by manipulation of the anisotropy in the Josephson coupling strength. 

We can define the fluxon mass in our ladder by equating the charging energy E c = 
C/2J2ij Vij t° the kinetic energy of a soliton of mass M v : Eki n = \M v v% 0. Since E c can 
be directly calculated in our simulation, while v v can be calculated from (V), this gives 
an unambiguous way to determine M v . For rj 2 = 0.5, we find E c /C ~ 110((\^) /I C R) 2 , 
in the flux- flow regime (region I of Fig. 1). This gives ~ 3.4C0q/o 2 , more than six 
times the usual estimate for the vortex mass in a 2D square lattice [13[]. Similarly, the 
vortex friction coefficient 7 can be estimated by equating the rate of energy dissipation, 
Edis = 1/2 Ynj Vij/Riji to This estimate yields 7 7 ~ 3A(f>Q/(Ra 2 ), once again more 

than six times the value predicted for 2D arrays [/J. This large dissipation explains the 
absence of ballistic motion for this anisotropy At larger values rjj = 5 and 2.5, a 

similar calculation gives AfJ ~ 0.28 and 0.34$/ {Ra 2 ), 7' ~ 0.28 and 0.34<f) 2 /{Ra 2 ). These 
lower values of 'j 1 , but especially the low pinning energies, may explain why ballistic motion 
is possible at these values of rjj. 

7 



To conclude, we have shown analytically that the dynamics of an underdamped Joseph- 
son ladder can be approximately reduced to those of a discrete sine-Gordon equation. As 
confirmation, we showed numerically that the ladder exhibits many phenomena previously 
reported in a discrete sine-Gordon system. These include separate flux-flow and whirling 
regimes in the IV characteristics, voltage steps in the flux flow regime which arise from lock- 
ing of vortices to phason excitations on the ladder, and ballistic vortex propagation on the 
steps at sufficiently high anisotropy. It would be of interest to construct anisotropic ladders 
in order to verify some of these predictions, and to attempt to extend these arguments to 
2D arrays. 
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supported by NSF Grant DMR94-02131 and by the Midwest Superconductivity Consortium 
through DOE Grant DE-FG02-90ER-45427. SR was supported by Ohio State University 
Postdoctoral Fellowship. 
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FIGURES 

FIG. 1. Calculated IV curves for various values of / in an 8 x 1 ladder ring. / = Ba 2 /4>q = 1/8 
is the number of flux quanta per plaquette. We use (3 = 33 and t]j = 0.71. Plotted voltages are 
divided by /. Inset: schematic of ring topology used in simulation. A uniform current is injected 
into the inner grains and drawn out from the outer grains. 

FIG. 2. Calculated IV curves for / = 1/8, (3 = 61, and several values of the anisotropy 
parameter rjj. The inset shows the voltages of steps corresponding to locking number n determined 
from our numerical result compared to those calculated from eq. 4 (solid lines). 

FIG. 3. Local phase difference ifi(t) (solid lines) and voltage Vi(t) (Broken lines) for i = 0, 
/ = 1/8, r]j = 1.25, (3 = 61, and three values of the applied current: I/I c = 0.22 and 0.41, 
corresponding to currents on voltage plateaus with n = 3, 2; and I/I c = 0.8, corresponding to a 
current in the whirling regime. Note the different scales in the axes for the whirling regime. 

FIG. 4. Plot of Vo(t) and dVo(t)/dt at x = after driving current is turned off at time t = 0. 
We use rjj = 5, and I/I c = 0.41, corresponding to an n = 1 voltage plateau. 
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